// 
// File:   stale.h
// Author: kajetan
//
// Created on 18 listopad 2007, 22:21
//

#ifndef _STALE_H
#define	_STALE_H

const double multiplier = 1.;

const double THETA  = 0.0;          // oddzialywanie z woda [?]

//const double D      = 3.4e-10;      // [m]
//const double EG     = -1.51;        // [eV]
//const double EA     = -0.95;        // [eV]
//const double ALPHAG = 1.587e-10;    // [m]
//const double ALPHAA = 2.e-10;       // [m]

const int G      = 0;
const int C      = 1;
const int A      = 2;
const int T      = 3;

// Hartree - eV 
const double Eh     = 27.2113845;

const double Em     = 2.71828183;

const double hbar      = 1.;

/*  
 * sprzezenia podane w eV
 */
const double sprz[][4] = 
{
    { 0.084 / Eh, 0.110 / Eh, 0.089 / Eh, 0.137 / Eh },
    { 0.042 / Eh, 0.041 / Eh, 0.029 / Eh, 0.100 / Eh },
    { 0.049 / Eh, 0.061 / Eh, 0.030 / Eh, 0.105 / Eh },
    { 0.085 / Eh, 0.076 / Eh, 0.086 / Eh, 0.158 / Eh }         
};  

const double E[4] = { -1.51 / Eh, -0.56 / Eh, -0.95 / Eh, -0.97 / Eh };

const double Pi     = 3.1415926535897;

const double K      = 8.61734e-5;   // stala Boltzmanna [eV/K]
const double Temp   = 300.;         // temperatura [K]

const double sigmaSq        = 4.;   // [j.a.]

const int KROKOW_NA_WYDRUK  = 1;
const int ELEMENTY          = 50;
const int KROKOW            = 1000000;

const double DT     = 10.;                       // timestep [j.a.]
const double DST    = 8.;                       // odl. miedzy zasadami [j.a.]
    
const double limitAkceptor  = (0.99);

// VACUUM
const double ampV[]  = { 0.15 / 0.52917721,      0.8 / 0.52917721};        // amplituda drgan [j.a.]
const double freqV[] = { 250.0e12 * 2.41888e-17, 1.5e12 * 2.41888e-17 };   // czestosc drgan [j.a.-1]

// WATER
const double ampW[]  = { 0.15 / 0.52917721,      0.8 / 0.52917721,       0.5 / 0.52917721};        // amplituda drgan [j.a.]
const double freqW[] = { 250.0e12 * 2.41888e-17, 2.0e12 * 2.41888e-17,   0.2e12 * 2.41888e-17};        // czestosc drgan [j.a.-1]

// const double amp    = 0.8 * 0.52917721;      // amplituda drgan [j.a.]
//const double freq   = 2.5e12;                 // czestosc drgan [s^-1]

/*
 *  1 j.a. czasu - 2.4188e-17 s (h/(2 pi E_{h})) - 2.4188e-5 ps
 *  
 */


#endif	/* _STALE_H */

